Enhanced RBE of Particle Radiation Depends on Beam Size in the Micrometer Range.

High-linear energy transfer (LET) radiation, such as heavy ions is associated with a higher relative biological effectiveness (RBE) than low-LET radiation, such as photons. Irradiation with low- and high-LET particles differ in the interaction with the cellular matter and therefore in the spatial dose distribution. When a single high-LET particle interacts with matter, it results in doses of up to thousands of gray (Gy) locally concentrated around the ion trajectory, whereas the mean dose averaged over the target, such as a cell nucleus is only in the range of a Gy. DNA damage therefore accumulates in this small volume. In contrast, up to hundreds of low-LET particle hits are required to achieve the same mean dose, resulting in a quasi-homogeneous damage distribution throughout the cell nucleus. In this study, we investigated the dependence of RBE from different spatial dose depositions using different focused beam spot sizes of proton radiation with respect to the induction of chromosome aberrations and clonogenic cell survival. Human-hamster hybrid (AL) as well as Chinese hamster ovary cells (CHO-K1) were irradiated with focused low LET protons of 20 MeV (LET = 2.6 keV/µm) beam energy with a mean dose of 1.7 Gy in a quadratic matrix pattern with point spacing of 5.4 × 5.4 µm2 and 117 protons per matrix point at the ion microbeam SNAKE using different beam spot sizes between 0.8 µm and 2.8 µm (full width at half maximum). The dose-response curves of X-ray reference radiation were used to determine the RBE after a 1.7 Gy dose of radiation. The RBE for the induction of dicentric chromosomes and cell inactivation was increased after irradiation with the smallest beam spot diameter (0.8 µm for chromosome aberration experiments and 1.0 µm for cell survival experiments) compared to homogeneous proton radiation but was still below the RBE of a corresponding high LET single ion hit. By increasing the spot size to 1.6-1.8 µm, the RBE decreased but was still higher than for homogeneously distributed protons. By further increasing the spot size to 2.7-2.8 µm, the RBE was no longer different from the homogeneous radiation. Our experiments demonstrate that varying spot size of low-LET radiation gradually modifies the RBE. This underlines that a substantial fraction of enhanced RBE originates from inhomogeneous energy concentrations on the µm scale (mean intertrack distances of low-LET particles below 0.1 µm) and quantifies the link between such energy concentration and RBE. The missing fraction of RBE enhancement when comparing with high-LET ions is attributed to the high inner track energy deposition on the nanometer scale. The results are compared with model results of PARTRAC and LEM for chromosomal aberration and cell survival, respectively, which suggest mechanistic interpretations of the observed radiation effects.

Intercomparison of dose enhancement ratio and secondary electron spectra for gold nanoparticles irradiated by X-rays calculated using multiple Monte Carlo simulation codes.

PURPOSE: Targeted radiation therapy has seen an increased interest in the past decade. In vitro and in vivo experiments showed enhanced radiation doses due to gold nanoparticles (GNPs) to tumors in mice and demonstrated a high potential for clinical application. However, finding a functionalized molecular formulation for actively targeting GNPs in tumor cells is challenging. Furthermore, the enhanced energy deposition by secondary electrons around GNPs, particularly by short-ranged Auger electrons is difficult to measure. Computational models, such as Monte Carlo (MC) radiation transport codes, have been used to estimate the physical quantities and effects of GNPs. However, as these codes differ from one to another, the reliability of physical and dosimetric quantities needs to be established at cellular and molecular levels, so that the subsequent biological effects can be assessed quantitatively. METHODS: In this work, irradiation of single GNPs of 50 nm and 100 nm diameter by X-ray spectra generated by 50 and 100 peak kilovoltages was simulated for a defined geometry setup, by applying multiple MC codes in the EURADOS framework. RESULTS: The mean dose enhancement ratio of the first 10 nm-thick water shell around a 100 nm GNP ranges from 400 for 100 kVp X-rays to 600 for 50 kVp X-rays with large uncertainty factors up to 2.3. CONCLUSIONS: It is concluded that the absolute dose enhancement effects have large uncertainties and need an inter-code intercomparison for a high quality assurance; relative properties may be a better measure until more experimental data is available to constrain the models.

MODELING STUDIES ON DICENTRICS INDUCTION AFTER SUB-MICROMETER FOCUSED ION BEAM GRID IRRADIATION.

The biophysical simulation tool PARTRAC contains modules for DNA damage response representing non-homologous end joining of DNA double-strand breaks (DSB) and the formation of chromosomal aberrations. Individual DNA ends from the induced DSB are followed regarding both their enzymatic processing and spatial mobility, as is needed for chromosome aberrations to arise via ligating broken ends from different chromosomes. In particular, by tracking the genomic locations of the ligated fragments and the positions of centromeres, the induction of dicentrics can be modeled. In recent experiments, the impact of spatial clustering of DNA damage on dicentric yields has been assessed in AL human-hamster hybrid cells: Defined numbers of 20 MeV protons (linear energy transfer, LET 2.6 keV/µm), 45 MeV Li ions (60 keV/µm) and 55 MeV C ions (310 keV/µm) focused to sub-µm spot sizes were applied with the ion microbeam SNAKE in diverse grid modes, keeping the absorbed dose constant. The impact of the µm-scaled spatial distribution of DSB (focusing effect) has thus been separated from nm-scaled DSB complexity (LET effect). The data provide a unique benchmark for the model calculations. Model and parameter refinements are described that enabled the simulations to largely reproduce both the LET-dependence and the focusing effect as well as the usual biphasic rejoining kinetics. The predictive power of the refined model has been benchmarked against dicentric yields for photon irradiation.

MODELING DNA DAMAGE BY PHOTONS AND LIGHT IONS OVER ENERGY RANGES USED IN MEDICAL APPLICATIONS.

Comprehensive track structure-based simulations of DNA damage induced in human cells by photons (5 keV-1.3 MeV) and light ions (0.25-512 MeV/u) were performed with PARTRAC. DNA strand breaks, double-strand breaks and their clustering were scored. Effective LET values were established for photons that provide LET-dependent damage yields in agreement with the data for ions. The resulting database captures the variations of biological effectiveness with radiation quality. In particular, it can help compare the effectiveness of conventional radiotherapy using photon beams with techniques relying on proton or ion beams.

WHAT ROLES FOR TRACK-STRUCTURE AND MICRODOSIMETRY IN THE ERA OF -omics AND SYSTEMS BIOLOGY?

Ionizing radiation is a peculiar perturbation when it comes to damage to biological systems: it proceeds through discrete energy depositions, over a short temporal scale and a spatial scale critical for subcellular targets as DNA, whose damage complexity determines the outcome of the exposure. This lies at the basis of the success of track structure (and nanodosimetry) and microdosimetry in radiation biology. However, such reductionist approaches cannot account for the complex network of interactions regulating the overall response of the system to radiation, particularly when effects are manifest at the supracellular level and involve long times. Systems radiation biology is increasingly gaining ground, but the gap between reductionist and holistic approaches is becoming larger. This paper presents considerations on what roles track structure and microdosimetry can have in the attempt to fill this gap, and on how they can be further exploited to interpret radiobiological data and inform systemic approaches.

A New Standard DNA Damage (SDD) Data Format.

Our understanding of radiation-induced cellular damage has greatly improved over the past few decades. Despite this progress, there are still many obstacles to fully understand how radiation interacts with biologically relevant cellular components, such as DNA, to cause observable end points such as cell killing. Damage in DNA is identified as a major route of cell killing. One hurdle when modeling biological effects is the difficulty in directly comparing results generated by members of different research groups. Multiple Monte Carlo codes have been developed to simulate damage induction at the DNA scale, while at the same time various groups have developed models that describe DNA repair processes with varying levels of detail. These repair models are intrinsically linked to the damage model employed in their development, making it difficult to disentangle systematic effects in either part of the modeling chain. These modeling chains typically consist of track-structure Monte Carlo simulations of the physical interactions creating direct damages to DNA, followed by simulations of the production and initial reactions of chemical species causing so-called "indirect" damages. After the induction of DNA damage, DNA repair models combine the simulated damage patterns with biological models to determine the biological consequences of the damage. To date, the effect of the environment, such as molecular oxygen (normoxic vs. hypoxic), has been poorly considered. We propose a new standard DNA damage (SDD) data format to unify the interface between the simulation of damage induction in DNA and the biological modeling of DNA repair processes, and introduce the effect of the environment (molecular oxygen or other compounds) as a flexible parameter. Such a standard greatly facilitates inter-model comparisons, providing an ideal environment to tease out model assumptions and identify persistent, underlying mechanisms. Through inter-model comparisons, this unified standard has the potential to greatly advance our understanding of the underlying mechanisms of radiation-induced DNA damage and the resulting observable biological effects when radiation parameters and/or environmental conditions change.

AT THE PHYSICS-BIOLOGY INTERFACE: THE NEUTRON AFFAIR.

We present predictions of neutron relative biological effectiveness (RBE) for cell irradiations with neutron beams at PTB-Braunschweig. A neutron RBE model is adopted to evaluate initial DNA damage induction given the neutron-induced charged particle field. RBE values are predicted for cell exposures to quasi-monoenergetic beams (0.56 MeV, 1.2 MeV) and to a broad energy distribution neutron field with dose-averaged energy of 5.75 MeV. Results are compared to what obtained with our RBE predictions for neutrons at similar energies, when a 30-cm sphere is irradiated in an isotropic neutron field. RBE values for experimental conditions are higher for the lowest neutron energies, because, as expected, target geometry determines the weight of the low-effectiveness photon component of the neutron dose. These results highlight the importance of characterizing neutron fields in terms of physical interactions, to fully understand neutron-induced biological effects, contributing to risk estimation and to the improvement of radiation protection standards.

Comprehensive track-structure based evaluation of DNA damage by light ions from radiotherapy-relevant energies down to stopping.

Track structures and resulting DNA damage in human cells have been simulated for hydrogen, helium, carbon, nitrogen, oxygen and neon ions with 0.25-256 MeV/u energy. The needed ion interaction cross sections have been scaled from those of hydrogen; Barkas scaling formula has been refined, extending its applicability down to about 10 keV/u, and validated against established stopping power data. Linear energy transfer (LET) has been scored from energy deposits in a cell nucleus; for very low-energy ions, it has been defined locally within thin slabs. The simulations show that protons and helium ions induce more DNA damage than heavier ions do at the same LET. With increasing LET, less DNA strand breaks are formed per unit dose, but due to their clustering the yields of double-strand breaks (DSB) increase, up to saturation around 300 keV/µm. Also individual DSB tend to cluster; DSB clusters peak around 500 keV/µm, while DSB multiplicities per cluster steadily increase with LET. Remarkably similar to patterns known from cell survival studies, LET-dependencies with pronounced maxima around 100-200 keV/µm occur on nanometre scale for sites that contain one or more DSB, and on micrometre scale for megabasepair-sized DNA fragments.

The origin of neutron biological effectiveness as a function of energy.

The understanding of the impact of radiation quality in early and late responses of biological targets to ionizing radiation exposure necessarily grounds on the results of mechanistic studies starting from physical interactions. This is particularly true when, already at the physical stage, the radiation field is mixed, as it is the case for neutron exposure. Neutron Relative Biological Effectiveness (RBE) is energy dependent, maximal for energies ~1 MeV, varying significantly among different experiments. The aim of this work is to shed light on neutron biological effectiveness as a function of field characteristics, with a comprehensive modeling approach: this brings together transport calculations of neutrons through matter (with the code PHITS) and the predictive power of the biophysical track structure code PARTRAC in terms of DNA damage evaluation. Two different energy dependent neutron RBE models are proposed: the first is phenomenological and based only on the characterization of linear energy transfer on a microscopic scale; the second is purely ab-initio and based on the induction of complex DNA damage. Results for the two models are compared and found in good qualitative agreement with current standards for radiation protection factors, which are agreed upon on the basis of RBE data.

Sub-micrometer 20MeV protons or 45MeV lithium spot irradiation enhances yields of dicentric chromosomes due to clustering of DNA double-strand breaks.

In conventional experiments on biological effects of radiation types of diverse quality, micrometer-scale double-strand break (DSB) clustering is inherently interlinked with clustering of energy deposition events on nanometer scale relevant for DSB induction. Due to this limitation, the role of the micrometer and nanometer scales in diverse biological endpoints cannot be fully separated. To address this issue, hybrid human-hamster AL cells have been irradiated with 45MeV (60keV/µm) lithium ions or 20MeV (2.6keV/µm) protons quasi-homogeneously distributed or focused to 0.5×1µm(2) spots on regular matrix patterns (point distances up to 10.6×10.6µm), with pre-defined particle numbers per spot to provide the same mean dose of 1.7Gy. The yields of dicentrics and their distribution among cells have been scored. In parallel, track-structure based simulations of DSB induction and chromosome aberration formation with PARTRAC have been performed. The results show that the sub-micrometer beam focusing does not enhance DSB yields, but significantly affects the DSB distribution within the nucleus and increases the chance to form DSB pairs in close proximity, which may lead to increased yields of chromosome aberrations. Indeed, the experiments show that focusing 20 lithium ions or 451 protons per spot on a 10.6µm grid induces two or three times more dicentrics, respectively, than a quasi-homogenous irradiation. The simulations reproduce the data in part, but in part suggest more complex behavior such as saturation or overkill not seen in the experiments. The direct experimental demonstration that sub-micrometer clustering of DSB plays a critical role in the induction of dicentrics improves the knowledge on the mechanisms by which these lethal lesions arise, and indicates how the assumptions of the biophysical model could be improved. It also provides a better understanding of the increased biological effectiveness of high-LET radiation.

Simulation on the molecular radiosensitization effect of gold nanoparticles in cells irradiated by x-rays.

Abundant studies have focused on the radiosensitization effect of gold nanoparticles (GNPs) in the cellular environment with x-ray irradiation. To better understand the physical foundation and to initially study the molecular radiosensitization effect within the nucleus, a simple cell model with detailed DNA structure in the central nucleus was set up and complemented with different distributions of single and multiple GNPs in this work. With the biophysical Monte Carlo simulation code PARTRAC, the radiosensitization effects on both physical quantities and primary biological responses (DNA strand breaks) were simulated. The ratios of results under situations with GNPs compared to those without GNPs were defined as the enhancement factors (EFs). The simulation results show that the presence of GNP can cause a notable enhancement effect on the energy deposition within a few micrometers from the border of GNP. The greatest upshot appears around the border and is mostly dominated by Auger electrons. The enhancement effect on the DNA strand breakage becomes smaller because of the DNA distribution inside the nucleus, and the corresponding EFs are between 1 and 1.5. In the present simulation, multiple GNPs on the nucleus surface, the 60 kVp x-ray spectrum and the diameter of 100 nm are relatively more effective conditions for both physical and biological radiosensitization effects. These results preliminarily indicate that GNP can be a good radiosensitizer in x-ray radiotherapy. Nevertheless, further biological responses (repair process, cell survival, etc) need to be studied to give more accurate evaluation and practical proposal on GNP's application in clinical treatment.

Energy dependence of the complexity of DNA damage induced by carbon ions.

To assess the complexity of DNA damage induced by carbon ions as a function of their energy and LET, 2-Gy irradiations by 100 keV u(-1)-400 MeV u(-1) carbon ions were investigated using the PARTRAC code. The total number of fragments and the yield of fragments of <30 bp were calculated. The authors found a particularly important contribution of DNA fragmentation in the range of <1 kbp for specific energies of <6 MeV u(-1). They also considered the effect of different specific energies with the same LET, i.e. before and after the Bragg peak. As a first step towards a full characterisation of secondary particle production from carbon ions interacting with tissue, a comparison between DNA-damage induction by primary carbon ions and alpha particles resulting from carbon break-up is presented, for specific energies of >1 MeV u(-1).

Cross-section scaling for track structure simulations of low-energy ions in liquid water.

Radiation damage by low-energy ions significantly contributes to the high biological efficiency of ion beams in distal Bragg peak regions as well as to the energy-dependent efficiency of neutron irradiation. To enable assessing biological effects of ions at energies <1 MeV u(-1) with track-structure based models, a Barkas-like scaling procedure is developed that provides ion cross sections in liquid water based on those for hydrogen ions. The resulting stopping power and range for carbon ions agree with the ICRU 73 database and other low-energy stopping power data. The method represents the basis for extending PARTRAC simulations of light ion track structures and biological effects down to the keV u(-1) range.

Modelling proton bunches focussed to submicrometre scales: low-LET radiation damage in high-LET-like spatial structure.

Microbeam experiments approximating high-LET tracks by bunches of lower-LET particles focussed to submicrometre scales (Schmid et al. 2012, Phys. Med. Biol. 57, 5889) provide an unprecedented benchmark for models of biological effects of radiation. PARTRAC track structure-based Monte Carlo simulations have verified that focussed 20 MeV proton bunches resemble the radial dose distributions of single 55 MeV carbon ions as used in the experiments. However, the predicted yields of double-strand break and short (<1 kbp) DNA fragments by focussed protons correspond to homogeneous proton irradiation and are much smaller than for carbon tracks. The calculated yields of dicentrics overestimate the effect of focussing but reproduce the fourfold difference between carbon ions and homogeneously distributed protons. The extent to which focussed low-LET particles approximate high-LET radiation is limited by the achievable focussing: submicrometre focussing of proton bunches cannot reproduce local nanometre clustering, i.e. DNA damage complexity characteristic of high-LET radiation.

Chromosome aberration model combining radiation tracks, chromatin structure, DSB repair and chromatin mobility.

The module that simulates the kinetics and yields of radiation-induced chromosome aberrations within the biophysical code PARTRAC is described. Radiation track structures simulated by Monte Carlo methods are overlapped with multi-scale models of DNA and chromatin to assess the resulting DNA damage. Spatial mobility of individual DNA ends from double-strand breaks is modelled simultaneously with their processing by the non-homologous end-joining enzymes. To score diverse types of chromosome aberrations, the joined ends are classified regarding their original chromosomal location, orientation and the involvement of centromeres. A comparison with experimental data on dicentrics induced by gamma and alpha particles shows that their relative dose dependence is predicted correctly, although the absolute yields are overestimated. The critical model assumptions on chromatin mobility and on the initial damage recognition and chromatin remodelling steps and their future refinements to solve this issue are discussed.

Mechanistic modelling of radiation-induced bystander effects.

A model of radiation-induced bystander effects is presented that explicitly takes into account the transient nature of bystander signal emission post-irradiation, signal lifetime and the non-linear cellular response to the signals. Data are analysed on mutagenesis induced in human lymphoblasts in medium transfer experiments, in which the signal build-up time, medium dilution and the duration of reporter cells' exposure to the medium were varied. The model implies that the cellular release of bystander signals decreases rather slowly, with a characteristic time of about a day, whereas the signal itself decays with a lifetime of about an hour.

Impact of intercellular induction of apoptosis on low-dose radiation carcinogenesis.

In vitro data indicate that selective removal of oncogenic transformed cells by apoptosis induced via signalling by neighbouring cells may represent an important anti-carcinogenic process. Mechanistic modelling supports this concept and predicts that the phenomenon can stop the growth of a transformed cell population, forming a dormant pre-neoplastic lesion, or even remove the transformed clone completely. Radiation has been shown to enhance the underpinning signalling and increase the extent and rate of apoptosis induction in precancerous cells. Implications for low-dose radiation carcinogenesis are discussed based on in vitro data and mechanistic modelling. The possibility is outlined for radiation to act in a pro-carcinogenic manner, i.e. to reduce rather than enhance the removal of transformed cells by apoptosis. The effects of radiation exposure during early or late carcinogenesis are discussed.

Modeling dose deposition and DNA damage due to low-energy ß(-) emitters.

One of the main issues of low-energy internal emitters concerns the very short ranges of the beta particles, versus the dimensions of the biological targets. Depending on the chemical form, the radionuclide may be more concentrated either in the cytoplasm or in the nucleus of the target cell. Consequently, since in most cases conventional dosimetry neglects this issue it may overestimate or underestimate the dose to the nucleus and hence the biological effects. To assess the magnitude of these deviations and to provide a realistic evaluation of the localized energy deposition by low-energy internal emitters, the biophysical track-structure code PARTRAC was used to calculate nuclear doses, DNA damage yields and fragmentation patterns for different localizations of radionuclides in human interphase fibroblasts. The nuclides considered in the simulations were tritium and nickel-63, which emit electrons with average energies of 5.7 (range in water of 0.42 µm) and 17 keV (range of 5 µm), respectively, covering both very short and medium ranges of beta-decay products. The simulation results showed that the largest deviations from the conventional dosimetry occur for inhomogeneously distributed short-range emitters. For uniformly distributed radionuclides selectively in the cytoplasm but excluded from the cell nucleus, the dose in the nucleus is 15% of the average dose in the cell in the case of tritium but 64% for nickel-63. Also, the numbers of double-strand breaks (DSBs) and the distributions of DNA fragments depend on subcellular localization of the radionuclides. In the low- and medium-dose regions investigated here, DSB numbers are proportional to the nuclear dose, with about 50 DSB/Gy for both studied nuclides. In addition, DSB numbers on specific chromosomes depend on the radionuclide localization in the cell as well, with chromosomes located more peripherally in the cell nucleus being more damaged by short-ranged emitters in cytoplasm compared with chromosomes located more centrally. These results illustrate the potential for over- or underestimating the risk associated with low-energy emitters, particularly for tritium intake, when their distribution at subcellular levels is not appropriately considered.

Integration of Monte Carlo simulations with PFGE experimental data yields constant RBE of 2.3 for DNA double-strand break induction by nitrogen ions between 125 and 225 keV/µm LET.

The number of small radiation-induced DNA fragments can be heavily underestimated when determined from measurements of DNA mass fractions by gel electrophoresis, leading to a consequent underestimation of the initial DNA damage induction. In this study we reanalyzed the experimental results for DNA fragmentation and DNA double-strand break (DSB) yields in human fibroblasts irradiated with γ rays and nitrogen ion beams with linear energy transfer (LET) equal to 80, 125, 175 and 225 keV/µm, originally measured by Höglund et al. (Radiat Res 155, 818-825, 2001 and Int J Radiat Biol 76, 539-547, 2000). In that study the authors converted the measured distributions of fragment masses into DNA fragment distributions using mid-range values of the measured fragment length intervals, in particular they assumed fragments with lengths in the interval of 0-48 kbp had the mid-range value of 24 kbp. However, our recent detailed simulations with the Monte Carlo code PARTRAC, while reasonably in agreement with the mass distributions, indicate significantly increased yields of very short fragments by high-LET radiation, so that the actual average fragment lengths, in the interval 0-48 kbp, 2.4 kbp for 225 keV/µm nitrogen ions were much shorter than the assumed mid-range value of 24 kbp. When the measured distributions of fragment masses are converted into fragment distributions using the average fragment lengths calculated by PARTRAC, significantly higher yields of DSB related to short fragments were obtained and resulted in a constant relative biological effectiveness (RBE) for DSB induction yield of 2.3 for nitrogen ions at 125-225 keV/µm LET. The previously reported downward trend of the RBE values over this LET range for DSB induction appears to be an artifact of an inadequate average fragment length in the smallest interval.